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Starting from a two-orbital model and based on the random phase approximation, spin excitations 
in the superconducting state of the newly discovered BiS2 superconductors with three possible 
pairing symmetries are studied theoretically. We show that spin response is uniquely determined 
by the pairing symmetry. Possible spin resonance excitations might occur for the d-wave symmetry 
at an incommensurate momentum about (0.77T, Q.7tt). For the p-wave symmetry the transverse spin 
excitation near (0, 0) is enhanced. For the s-wave pairing symmetry there is no spin resonance 
signature. These distinct features may be used for probing or determining the pairing symmetry in 
this newly discovered compound. 
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Very recently, new BiS2 based superconducting (SC) 
materials were discovered 1-0]. The materials reported 
so far include B14O4S3 LaOi_ K F x BiS2 [2j, and 
NdOi_ a; F :I ;BiS2 3]. These materials share some com- 
mon features with cuprate-based and iron-based high-Tc 
materials, e.g., they are layered materials with BiS2 lay- 
ers being the SC plane. Superconductivity is realized by 
adding defects into the parent systems. Thus this family 
might open a new door in the studies of unconventional 
superconductivity, providing a new insight into the origin 
of high-Tc superconductivity. Actually, it has already at- 
tracted considerable interest in investigating the physical 
properties of them ■ 

Understanding the pairing mechanism and probing the 
pairing symmetry are two of most important issues in 
the studies of SC materials. Now it is believed that the 
spin fluctuation should be responsible for the SC pair- 
ing in many families of unconventional SC materials, in- 
cluding heavy fermion materials, cuprates and iron-based 
ones [9(. For BiS2 based materials, in was found from 
band calculation that there exists a good nesting of Fermi 
surface, thus the spin fluctuation should also be consid- 
erably strong On the other hand, it was proposed 
in Ref. Q that the SC pairing is strong and exceeds the 
limit of the phonon mediated picture. Therefore, the 
spin fluctuation is a good candidate to be responsible for 
the SC pairing in this family. In the mean time, the 
results of the spin susceptibility are usually phase sen- 
sitive due to the existence of coherence factor, thus it 
has been a p owerful tool for probing the pairing symme- 
tries [ig - ll3j . At last, many exotic features, e.g., the spin 
resonance and possible incommensurate (IC) spin excita- 
tion, have attracted broad interest in previous studies of 
various families of unconventional superconductors [hH - 
[l6j . Therefore, at this stage it is highly desirable and 
of importance to study the spin excitations in the BiS2- 
based materials and to look into the pairing symmetry 
theoretically. 

Band structure for LaOi_ 3; F a ;BiS2 sample was ob- 



tained from first principles calculation in Ref. Q ■ A sim- 
plified two-orbital (considering Bi 6p orbital) model was 
proposed by neglecting the interlayer coupling and the 
bands far away from Fermi surface In this paper, by 
employing this two-orbital model, we study theoretically 
the spin fluctuations based on the random phase approx- 
imation (RPA). Three possible SC pairing symmetries 
(e?-wave, extended s-wave, and p-wave) have been con- 
sidered to study the spin response in the SC state. We 
find that the spin response is enhanced at an IC mo- 
mentum near (71", 7r) [or (0,0)] for the <i-wave [or p-wave] 
pairing symmetry, while there is no significant feature of 
the spin excitations for the s-wave pairing. These could 
be used to test the pairing symmetry experimentally. 
Our starting model is expressed as, 

where /i, v are the orbital indices. In the present work, 
we use the two orbital model with the hopping constants 
being from Ref. Q. Ha is the SC pairing term. The 
last term is the on-site electron interaction and will take 
effect in the RPA framework. Usually the interaction 
term should include the onsite repulsive interaction, and 
the Hund coupling term. Both terms play a similar role 
and enhance the spin excitation after the RPA correc- 
tion. Thus we here consider only the onsite intraorbital 
interaction term. We have also checked numerically that 
the main results do not change when an additional inter- 
orbital interaction is added. 

We define the single-particle normal Green's 
function G and anomalous Green's function F(F) 
with G CT (k, ir) M „ = -(T T c M(T (k, -ir)4 (T (k, 0)), 



F ff (k,ir)^ 



-(T T c M<7 (k, -ir)c !/) _ CT (— k, 0)), and 



F°(k,zr) uu = -(T T cl a (k,-iT)c^_ a (-k,0)). 

Then the bare spin susceptibility can be defined as, 

= - w T,kA G t 3 ( k ) G fj( k + <?) 

-vo'fe(k)FZ'{k + q)], (2) 
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FIG. 1: (Color online) (a) The normal state Fermi surface, (b) 
The normal state spin excitations along different momentum 
cut. 



Here k = (k,iw„) and q = (q,iu). A(k) (A = G,F,F) 
are the Fourier transformations for A(k, it), a' = a for 
Xg 2 and a 1 = —a for Xo ' respectively. 

The anomalous Green's function in Eq.(2) satisfy 
Ffj(k) = F^i-k) = -F[i(k). For the spin singlet 
pairing, F^(k) — —F^ a (k), thus the formulas for \ zz 
and x^ are the same. For the triplet pairing, we have 
F^j(k) = F^ a (k), so the spin excitations for \ zz an d x H 
are different with an additional minus sign in front of the 
anomalous part of spin susceptibility. 

The spin susceptibility within RPA is given by 



xo(q,^) 



I - c/ixo(q,w) 



(3) 



where I is a unit matrix. The total (physical) spin sus- 
ceptibility x is YsijXir 

Before we present our numerical results, we need to 
determine the doping density. For the B14O4S3 sample, 
the band calculations suggest that the electron doping 
is x — 0.5 li! . For LaOi- 2; F l; BiS2 sample, supercon- 
ductivity is observed for x > 0.5 and it was proposed 
that the optimal doping density is slightly above 0.5 
For the NdOi_ a; F :c BiS2 sample, the range of SC state 
seems to be broad with superconductivity being observed 
for 0.1 < x < 0.7 [3|. We also noted that the Fermi 
surface topology changes significantly as the doping in- 
creases from 0.25 to 0.5 [J], thus the pairing symmetry 
and the spin excitations might be quite different for dif- 
ferent doping densities for the NdOi_ x F x BiS2 sample, 
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FIG. 2: (Color online) Three possible SC pairing, (a) Ak oc 
cosfca;— cosky. (b) Ak oc cos fc^+cos k y . (c) Ak oc sin(k x ^k y ). 



while this issue was not concerned in the present work. 
We focus our studies on the optimal doped sample with 
x = 0.55 used throughout the present work. Other pa- 
rameters are chosen as U — 1.5, Ao = 0.06. 

The normal state Fermi surface and spin susceptibil- 
ity as a function of momentum q are plotted in Fig. 1(a) 
and Fig. 1(b). As is seen, the spin excitations are mainly 
concentrated near the momentum (0,0) and {tt,tt). We 
also find that the low energy spin excitations are incom- 
mensurate, with the maximum spin excitations being at 
(0, ±6) and (it ± S, n ± S). The origin of these features is 
explained well based on the Fermi surface nesting picture. 
As seen in Fig. 1(a), since the band structure includes two 
quasi-one dimensional band, the Fermi surface exhibits a 
good nesting feature with the vectors shown in Fig. 1(a) 
by arrows [4(. 

The pairing symmetry can be estimated by assum- 
ing that SC pairing is mediated by the spin fluctua- 
tion [17T - [l9| . Generally, if the pairing potential around 
the momentum q determine the SC pairing, the rela- 
tion A(k) = — A(k + q) should be satisfied with k 
being the wave vector on the Fermi surface. With a 
weak-coupling analysis fl§| , we can obtain three possi- 
ble pairing symmetries, displayed in Fig. 2. For a sin- 
glet pairing symmetry, their exist two possible pair- 
ing symmetries, as shown in Fig. 2(a) and 2(b). The 
rf-wave symmetry with A(k) = Ao(cosfca; — cosfc y )/2 
shown in Fig. 2(a) is caused by the (7r,7r)-spin excitation. 
Meanwhile, both (0, ±5) and (ir, 7r)-spin excitation would 
give rise to the extended s-wave pairing symmetry with 
A(k) = Ao^osfcr+cosfcy)^. For the triplet pairing, the 
p-wave pairing with Ak = Aosin^ =p k y ) (with ' — ,+' 
for orbital 1,2, respectively) is favored. 

The numerical results of the spin excitations for the 
d-w&ve pairing symmetry are presented in Fig. 3. At 
the (7r, 7r) momentum the spin susceptibility increases as 
the frequency increases and no peak structure is seen. 
As the momentum is transferred to the IC one, i.e., 
q = (tt±5, tt±S) (5 is the incommensurability), the peak 
appears at the frequency uj = uj p . Our numerical results 
also show that u p increases as the incommensurability 5 
increases (not shown here). In the meantime, the inten- 
sities of the IC peaks also increase as 8 increases. At the 
momentum about (0.77T, 0.77r), the IC peaks reach the 
maximum intensity at the frequency about oj p = 0.08. 
At higher frequencies, spin excitations at the IC momen- 
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FIG. 3: (Color online) The imaginary parts of the spin sus- 
ceptibility as a function of frequency at different momenta for 
the d-wave SC state. The normal state spin excitations are 
also displayed for comparison. The inset of panel (a) shows 
the bare spin susceptibility. 



turn decreases as the frequency increases. As u > 0.1, 
the commensurate (ir, tt) spin excitation is dominant over 
IC ones, as seen in Fig. 3(a). In the meantime, we also 
looked into the spin excitations around the parallel direc- 
tion with q = (tt, 7t ± 8), and found that the parallel spin 
excitations are always smaller than the diagonal ones, 
therefore they are not presented here. 

We also plotted the RPA normal state spin susceptibil- 
ity at the momentum (0.77T, 0.7n) in Fig. 3(a). Obviously, 
the spin excitation is enhanced significantly in the SC 
state, such that it is quite like a resonance behavior. Usu- 
ally the spin resonance occurs for the pole condition of 
PRA factor and negligibly small bare spin susceptibility. 
To check whether it is the resonance or quasi-resonance 
behavior, we plotted the bare spin susceptibility in the 
inset of Fig. 3(a). As is seen, at very low energies the 
imaginary part of bare spin susceptibility is zero due to 
the presence of the spin gap. At the gap edge, there is 
a steplike rise of Im^o- Due to the Kramers-Kroenig re- 
lations the Rexo develops a peak structure. As a result, 
the spin resonance would be expected when 1 — UKe\o 
equals to zero (or small, corresponding to a small onsite 
interaction U). Actually for the present case we can con- 
clude that it is indeed a resonance behavior while U is 
not very large here so that the resonance peak is not very 
strong. It will increase greatly if a larger U is taken into 
account, but it is not presented in the present work. 

For small momenta, the spin excitations are incom- 
mensurate for all of the energies {uj < 0.3) we considered, 



(b) 




CO 

FIG. 4: (Color online) The imaginary parts of the spin sus- 
ceptibility as a function of frequency at different momenta for 
the s-wave SC state. 



with the IC peaks along the k y — direction. We plot- 
ted the frequency dependence of RPA spin susceptibility 
Imx(q, cj) for the momentum q = (0.067T, 0) in Fig. 3(b). 
As is seen, the intensity of the IC peak is suppressed sig- 
nificantly by the SC order. This is due to the presence of 
the coherence factor in the bare spin susceptibility fllT - 
[l3j]. In principle, for spin-singlet pairing symmetry, the 
spin response is enhanced if AkAk+ q < is satisfied, 
otherwise it would be suppressed in the SC state. Thus 
it is quite reasonable that the spin response is enhanced 
near (tt, tt) and suppressed around small momenta for the 
(i-wave pairing symmetry. 

The spin excitations for the extended s-wave are pre- 
sented in Fig. 4. Fig. 4(a) shows the spin excitation near 
(7T, tt). As is shown, the spin susceptibility increases al- 
most linearly as the frequencies at the commensurate mo- 
mentum. For the IC momentum (0.7tt, 0.77r), the spin 
susceptibility exhibits a broad peak. However, at this 
momentum the peak intensity nearly equals to that of 
the normal state shown in Fig. 3(a). Thus there is no 
spin resonance excitations. We also examine other mo- 
menta along various directions and no obvious resonant 
spin excitation was found. On the other hand, the IC 
behavior is similar to the case of d-wave symmetry. The 
spin response is incommensurate for low frequencies. The 
incommensurability decreases as the frequency increases 
and the spin excitations shift to the commensurate one 
at the frequency about 0.15. 

The spin excitation around a small momentum is dis- 
played in Fig. 4(b). The spin excitations do not show 
much difference compared to the normal state one. Ac- 
tually, since the nodal lines are close to the Fermi surface, 
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FIG. 5: (Color online) The imaginary parts of the spin sus- 
ceptibility as a function of frequency at different momenta for 
the p-wave SC state. The inset of panel (b) shows the bare 
spin susceptibility. 



the pairing gap at the fermi surface is quite small. As a 
result, there is nearly no significant features for the SC 
spin response. 

The spin excitations for the p-wave symmetry are de- 
picted in Fig. 5. In this case the transverse (x H ) and 
longitude (x zz ) spin susceptibility differ with the coher- 
ence factors. For both q = (tt,tt) and q = (0,0), one 
should have A k A k+q > 0. Thus generally x + ~ is greater 
than x zz ■ The enhancement of spin response in the SC 
state can only occur in the x^ channel. This can be 
seen clearly from Fig. 5 and the spin susceptibility for 
zz channel is suppressed by the SC pairing and always 
smaller than that for H — channel. Around the (tt, tt) mo- 
mentum, the spin excitation for x~* is similar to that of 



s-wave case, i.e., at low frequencies and the momentum 
q = (tt, tt), the spin susceptibility increases as frequency 
increases. Away from the commensurate momentum, a 
broad IC peak shows up. At higher frequencies, the com- 
mensurate spin excitation is dominant. 

Similar to the cases of d-wave and s-wave, the spin 
responses at small momenta are incommensurate for all 
of the frequencies. The strongest excitation occurs at 
the momentum (0.067r,0) (and its symmetric momenta). 
The spin excitation at this momentum is presented in 
Fig. 5(b). As is seen, the spin excitation in H — channel 
displays a peak structure at the frequency uj p = 0.058. 
The IC peak intensity is greater than that of normal state 
while it is not significant enough. Such enhanced spin ex- 
citation at this momentum can be explained through the 
bare spin excitation shown in the inset of Fig. 5(b). The 
origin of the enhancement is similar to that of spin reso- 
nance, with a step-like rise for Imxo and a peak structure 
for Rexo at the frequency lu p . The difference between the 
present case with the resonant spin excitation shown in 
Fig. I (a) is that here Im^o is finite at low frequencies, and 
thus the pole condition for the complex RPA factor was 
not satisfied. Therefore, it is a relatively weak feature, 
while we propose that if it is detected experimentally, 
it might be used as a signature to support the p-wave 
pairing symmetry. 

In summary, we have studied the spin excitations in 
the newly discovered BiS2 based superconductors based 
on a two orbital model. Assuming that the SC pairing 
is mediated by the spin fluctuation, three possible pair- 
ing symmetries are suggested. The spin excitations in 
the SC state is studied in detail. There exist some com- 
mon features for different pairing symmetries, such as 
the IC spin excitations at low frequencies. Meanwhile, 
distinct spin responses corresponding to different pairing 
symmetries can be used to probe the pairing symmetry 
experimentally. 
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